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Abstract 

We consider a model of socially interacting individuals that make a binary choice in a 
context of positive additive endogenous externalities. It encompasses as particular cases 
several models from the sociology and economics literature. We extend previous results 
to the case of a general distribution of idiosyncratic preferences, called here Idiosyncratic 
Willingnesses to Pay (IWP). 

Positive additive externalities yield a family of inverse demand curves that include the 
classical downward sloping ones but also new ones with non constant convexity. When j, 
the ratio of the social influence strength to the standard deviation of the IWP distribution, 
is small enough, the inverse demand is a classical monotonic (decreasing) function of the 
adoption rate. Even if the IWP distribution is mono-modal, there is a critical value of j 
above which the inverse demand is non monotonic, decreasing for small and high adoption 
rates, but increasing within some intermediate range. Depending on the price there are 
thus either one or two equilibria. 

Beyond this first result, we exhibit the generic properties of the boundaries limiting 
the regions where the system presents different types of equilibria (unique or multiple). 
These properties are shown to depend only on qualitative features of the IWP distribution: 
modality (number of maxima), smoothness and type of support (compact or infinite). The 
main results are summarized as phase diagrams in the space of the model parameters, on 
which the regions of multiple equilibria are precisely delimited. 
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1 Introduction 



1.1 Modeling social influences with heterogeneous agents 

There are many circumstances in social and economic contexts where, faced with different 
alternatives, the best choice for an individual depends on the choices of other individuals in the 
population. The decision of leaving a neighborhood [75J, to attend a seminar [75J or a crowded 
bar [2J to participate to collective actions such as strikes and riots [JB], are particular 
examples taken from social sciences. It has been suggested that social interactions may explain 
the school dropout [21] , the persistence in the educational level within some neighborhoods [30 
and the related consequences in the stratification of investment in human capital and economic 
segregation the large dispersion in urban crime through cities with similar characteristics 
|42j . the emergence of social norms |71j . the labor market behavior and related unemployment 
patterns [HH1 [EE], the housing demand [SB]) the existence of poverty traps [34 , the smoking 
behavior [S5J [55] [S3] , etc. 

Similarly, there is a growing economic literature that recognizes the influence on consumers 
of the social world they live in. In market situations like the subscription to a telephone 
network [5J [75J [301 123 or the choice of a computer operating system 52J , the willingness to 
pay generally depends not only on the individual preferences but also on the choice made by 
others [5TJ [75] . If the externality is positive the utility of the most popular choice increases even 
for individuals who otherwise would never make this choice. In other words, the conformity 
effect may dominate the heterogeneity of preferences, as pointed out by Bernheim [12J. General 
aspects of these issues have been discussed in the literature [10l[60]. Particular insightful papers 
are Becker's note [S] about restaurants pricing, and the qualitative analysis by Granoveter and 
Soong [17] of the consequences of interpersonal influences ("bandwagon effects" [571 I76| ) on the 
consumers demand and on the supply prices. 

In the present paper we consider the general properties of a model of socially interacting 
heterogeneous individuals that make a binary choice in a context of positive endogenous ex- 
ternalities. The model encompasses, as particular cases, most of the above mentioned models 
presented in the sociology and economics literature. In a forthcoming paper 44J we explore the 
consequences of the externalities on the economy, taking as an example the simplest market, 
i.e. that of a monopolist pricing a single good. 

In social sciences, the question of discrete (typically binary) choices with heterogeneous 
agents and positive externalities has been first addressed in the 70 's by Schelling [77] [75], who 
borrowed from Physics the concept of critical mass: in a repeated-decisions setting, depending 
on whether this critical mass is or not reached, the system may end up at very different equi- 
libria. Granovetter further develops Schelling's model, applying it to particular problems such 
as joining or not a riot [46], voting, etc [47] . The same topic is reconsidered within a statis- 
tical physics point of view in the early 80's by Galam et al [39]. The notion of critical mass 
is then related to the Physics concept of phase transition at a critical point, in the neighbor- 
hood of which the system may be extremely susceptible: by tipping effects, small microscopic 
changes can lead to drastic changes at the macroscopic level. Similar tools have been applied 
in 1980 by Kindermann and Snell |53j to the study of social networks. These authors intro- 
duced into the sociology and economics literature the equivalence between statistical physics 
approaches — that use the Boltzmann-Gibbs distribution — and Markov Random Fields. An- 
other physically-inspired approach for modeling social phenomena such as opinion diffusion 
has been developed by Weidlich and Haag [S31 IH2] in 1983, through a master equation and 
the Fokker-Planck approximation. Later, these physically inspired models of opinion contagion 
have been exploited in economics by Orlean [69 [ 170] for the analysis of mimetic behaviors in 
the context of financial markets. There is now a large and growing literature on opinion and 
innovation diffusion (see e.g. [5^1 EH1 [2E1 ED 1EI] ) closely related to the general discrete choice 
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model considered in this paper. Since the beginning of the 90's the general framework of social 
interactions in non-market contexts is reconsidered in a Beckerian way [5] 110] , in particular by 
Glaeser et al jUJ |42l SO] . 

The first application of statistical mechanics approaches in economics may be traced back 
to the pioneering work of Former [36] . Introducing an economic interpretation of the Ising 
model of ferromagnetism at finite temperature, he shows that strong externalities may hinder 
the stabilization of an economy. These models introduce Markov random fields (equivalently 
Boltzmann-Gibbs distributions) to model uncertainty in the decision making process, allowing 
for the definition of a general equilibrium concept. According to Former, Hildebrand's [48 justi- 
fication of the representative agent approach breaks down when agents' decisions are correlated 
due to their social interactions (for a discussion, see also [53]). 

A renewal of interest for models of binary decisions with externalities arose in economics 
in the 90's. On one side, Durlauf and collaborators [27j HH |29j [30] and Kirman and Weis- 
buch [94] among others, consider agents that choose an action according to a Boltzmann-Gibbs 
distribution, that is a logistic choice function, reflecting some random aspects in the agent's 
utility. In this context Brock [TU] and Blume [T21 H3] explicit the links between Game Theory 
and Statistical Mechanics, while Kirman and coworkers [57] show that the logistic choice func- 
tion may be seen as resulting from an exploration-exploitation compromise. These and other 
recent papers [H EU E21 [20l EU [40l EH EH [45j EH [16] analyze with statistical physics tools the 
consequences of positive social (market and non-market) interactions in the aggregate behavior 
of large populations (for a short introduction to statistical physics approaches see [33] and for 
their application to economics see [72]; see also [7] for a survey). Most of the above mentioned 
authors restrict the analysis of the model to the case where all the individuals have the same 
idiosyncratic preference. Heterogeneity in the population is introduced through the probabilis- 
tic decision-making process (random utility model [58, 59 ), like in Former's work |36j . Then, 
the actual equilibrium reached by the system depends on the fixed points of the decision dy- 
namics, generally a myopic best reply. An interesting characteristic of these models is that they 
present multiple equilibria for some range of the parameters. Becker [9] pointed out important 
consequences of these multiple equilibria, induced by externalities, on the economy: he suggests 
that they could be the reason of seemingly suboptimal pricing in situations of persistent excess 
demand. 

In this paper we consider intrinsically heterogeneous agents with fixed utilities, like in Mc- 
Fadden's approach to Quantal Choice models 62, 63 . Each individual has an Idiosyncratic 
Willingness to Pay (IWP) that remains fixed on time. We mainly (but not exclusively) decline 
the model within a market context, in which the binary choice corresponds to buying or not a 
given good at a posted price. This general setting allows us to generalize Becker's qualitative 
analysis [S] of the optimal pricing problem. Putting the price to zero allows us to recover the 
social sciences models. We assume that these IWP are random variables that are distributed 
among the population according to a given probability density function (pdf). 

We determine the possible equilibria of the system without assuming any precise decision- 
making dynamics. We show that the model's properties depend on the strength of the ex- 
ternality and on qualitative properties of the IWP pdf, like its modality class (the number 
of maxima), its smoothness properties and the kind of support. We display the main results 
on a plane whose axes arc the parameters of the model, namely, the average IWP and the 
strength of the social component, both measured in units of the standard deviation of the IWP 
distribution. Particular cases of our model have been published elsewhere 66, 45] . This paper 
extends those result to the case of a general IWP distribution. The particularly important case 
of a uni-modal pdf (with a single maximum) is thoroughly studied, but we also discuss the 
consequences of multi-modality. Our results are summarized on phase diagrams, that is, we 
plot in the parameters space the boundaries of the regions where different types of solutions 
exist . 
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Before entering into the details of our approach (section [2| we discuss the analogies and 
differences between our model and other models of social interactions (section 1.2 1 and we 
summarize our main results (section 1.3). 



1.2 More on related models 

In this section we briefly discuss the relationship between the model to be considered here 
and other models studied in the literature. Let us first consider models of discrete choices 
in the absence of externalities. According to the typology proposed by Anderson et al. [I], 
within the general framework of Random Utility Models (RUM) (SSJ [5j5] with additive stochastic 
utilities, there are two distinct approaches to individual choices: a "psychological" one and an 
"economic" one. In the psychological perspective (Thurstone [87], Luce [58]) the randomness is a 
time-dependent i.i.d. random variable: the random components of the idiosyncratic preferences 
are assumed to be independently drawn afresh by each individual from a given pdf, each time 
the choice has to be made. They are interpreted as individual temporary changes, or mistakes 
in the estimated utilities. In the simplest case — actually, the only one treated in the social 
and economic literature — the agents IWPs have identical deterministic parts and only differ 
by this random time- varying term which is systematically assumed to be drawn from a logistic 
pdf. In practice many approaches like in [20 consider the choice rule as deriving from a random 
utility model [55]. As shown by McFadden [55], in this context the logistic form is obtained if 
the random terms in the underlying Thurstone's discriminant process are i.i.d. Weibull random 
variables, i.e. have a double exponential (extreme value, type I) distribution (see also [1]). 

In the presence of strategic complementarities ( 22 ,23 ), the resulting model is well known 
in statistical physics: it corresponds to the standard Ising model, i.e. with ferromagnetic in- 
teractions and annealed disorder, that is, at finite temperature T . The latter is the inverse of 
the standard logistic parameter (3=1 /T and is thus proportional to the standard deviation of 
the IWP distribution. The ferromagnetic interaction constant J corresponds to the strength 
of the social externality. Introduced by Folmer [36 in the economics context, this standard 
Ising model has recently been reconsidered in the social and economic literature mainly by 
Durlauf and coworkers [3T] 1321 (US HI] and by Kirman and Weisbuch [MJ [57] . The correspond- 
ing equilibria are reminiscent of the Quantal Response Equilibria [64 used in the context of 
experimental economics and behavioral game theory. These are equilibria "on the average" , 
in the statistical sense (as in Physics): they do not correspond to the strict maximization of 
the utilities (that are random variables) but to that of their estimated or expected values. In 
the generally considered infinite population limit (where the variance of the expected values 
vanishes) the expected utilities are systematically smaller that the maximal ones. 

The standard Ising model is quite well understood |85| . Ising himself |51j gave an analytic 
description of its properties in the simplest case of a linear chain where each agent interacts only 
with his right and left nearest neighbors. There is also an analytical description of the stationary 
states of the model, due to Onsager [55] , in the case where the agents are situated on the vertices 
of a 2-dimcnsional square lattice, each having four neighbors. However, no analytic results exist 
for arbitrary neighborhoods except for the specific case of a global neighborhood, known as the 
mean field approximation in Physics. Accordingly, considering global neighborhoods, Brock and 
Durlauf [20] analyze the properties of the expected demand in the case of rational expectations 
under the assumption of a logistic distribution of such expectations (assuming thus double 
exponential random utilities). They find, in agreement with standard results in statistical 
mechanics |85j . that there exist either one, two or three solutions for the demand function, 
depending on the relative magnitudes of the idiosyncratic uniform social term, the variance of 
the stochastic term and the strength of the social effects. 

In the following we adopt instead McFadden's [53] economic approach (see also [5J0 H]): 
we assume that each agent has a willingness to pay invariable in time, that is different from 
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one agent to the other. In statistical physics this heterogeneity is called quenched disorder. 
The particular model we study is analogous to the ferromagnetic Random Field Ising Model 
(RFIM) at zero temperature (corresponding to the fact that the agents make deterministic 
choices). Thus, our modeling approach assumes the so called "risky" situation: an external 
observer (e.g. a seller) does not have access to the individual preferences, but may know their 
probability distribution. According to McFadden, "Thurstone's construction is appealing to an 
economist because the assumption that a single subject will draw independent utility functions 
in repeated choice settings and then proceed to maximize them is formally equivalent to a model 
in which the experimenter draws individuals randomly from a population with differing, but 
fixed, utility functions, and offers each a single choice; the latter model is consistent with the 
classical postulates of economic rationality" ([63], p 365). However, in the presence of social 
interactions this statement is in general incorrect. In a repeated choice setting, individual 
utilities evolve in time according to the others' decisions. In fact, the equilibria of systems 
reached through a dynamics that corresponds to an iterated game where agents make myopic 
choices at each time-step is one of the main concerns of statistical physics. It is known that 
the equilibria of systems with annealed or quenched disorder are of very different nature. The 
time average on a single agent and the population average at a given time do not necessarily 
coincide. 

In contrast to the standard Ising model at finite temperature (annealed disorder), the prop- 
erties of the RFIM with externalities, both at zero and at finite temperatures, are far from being 
fully understood. The properties of quenched disordered systems have been and are still the 
subject of numerous studies. Since the first studies of the RFIM, which date back to Aharony 
and Galam |38 , 37 , a number of important results have been published in the physics literature 
on this model (see e. g. [80]). Several variants of the RFIM have already been used in the 
context of socio-economic modeling, both by physicists and economists [39 ] l70 [ IT8 | l95j. 

The quenched-utilities model (RFIM at zero temperature) and the annealed-disorder model 
(standard mean field Ising model at finite temperature) have the same aggregate behavior 
(i.e. demand function for the market case) and equilibria under the following conditions - but 
essentially only under such conditions: 

1. the choice function with annealed utilities is identical to the cumulative distribution of 
the quenched IWPs; 

2. in the annealed case, equilibrium is reached through repeated best reply choices, where 
the expected demand is myopically estimated; 

3. the population size is infinite, guaranteeing that the variance of the demand vanishes in 
both models. 

However, the economic interpretation of these equilibria are very different: in the case of 
quenched utilities these are standard Nash equilibria, while in the case of annealed utilities 
these are similar (although not identical) to Quantal Response equilibria. 

1.3 Main specific results 

In the present paper we determine the equilibrium properties in the case of a global neighborhood 
with time invariant (quenched) random utilities, in the limit of an infinite number of agents. 
Since our paper focuses on equilibrium (static) properties, the social influence depends on the 
actual choices of the neighbors, in contrast with |20j . where the social influence in the surplus 
function depends on the agent's expected demand. 

Previous studies using annealed or quenched utilities consider specific probability distribu- 
tions, mostly a logistic or a Gaussian [87] . Some papers have determined conditions on the 
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choice function for having multi-equilibria [57J HO]- From the Physics literature we expect that 



specific properties near a critical point (a bifurcation point, see Section 3.2.1 below) are inde- 
pendent of the details of the model: this is used in [80] for describing the hysteresis effects in 
a family of (physical) systems at such a critical point, and exploited in 65] for the analysis of 
empirical socio-economic data in cases where the actual pdfs are not known. However, the full 
description of the phase diagram for an arbitrary pdf has not been done yet. Here we present 
this detailed analysis for a typical probability distribution of the IWP. We show how uniqueness 
or multiplicity of equilibria, related to convexity properties of the inverse demand functions, 
result from modality and smoothness of the pdf, as well as from the strength of the externality. 

We show that for small enough strengths of the social influence (the case of moderate social 
influence in [10]), the demand has a classical shape, that is, with a continuous decreasing 
adoption rate for increasing prices. However, if the ratio between the social influence strength 
and the standard deviation of the IWP distribution exceeds a critical value, the inverse demand 
function exhibits a non-classical, non-monotonic, behavior. As a result, depending on the price, 
there are either one or two stable equilibria for the demand: the positive (additive) externalities 
in a market context may give raise to a family of non-monotonic demand curves generalizing 
thus the classical ones. 

Beyond this first main result, we exhibit the generic properties of the boundaries limiting 
the regions where the system presents different types of equilibria (unique or multiple). We 
call these properties generic since we show that they depend only on qualitative features of 
the IWP distribution: modality (number of maxima), smoothness (continuity and derivability 
properties) and type of support (compact or infinite). The main results are summarized as 
phase diagrams in the space of the relevant parameters of the model, namely (i) the social 
influence strength and (ii) the difference between the population average of the IWP and the 
posted price, both parameters conveniently normalized by the standard deviation of the IWP 
distribution (which measures the importance of the heterogeneity). 

1.4 Organization of the paper 



The paper is organized as follows. In Sectio n [2| w e present the model: we first (section 2.1) 



specify the agents (customers) model, then in 2.2 we introduce a normalized form of the basic 



equations which is convenient for analyzing the demand, and in 2.3 we show on two simple 
extreme cases what to expect from these equations. In |2.4| we detail the families of probability 
distributions of the IWP covered by this paper. 

In Section [3] we analyze the aggregate demand (its collective behavior) for a generic smooth 
pdf. In |3 - 1 1 we introduce and study the direct and inverse demand functions; the demand 



phase diagram is derived in Section |3.2| in |3.2.1 we obtain the domain of multiple solutions 



which allows to plot the phase diagram (section 3.2.3). We analyze in details the vicinity of 
the bifurcation point (section |3.2.4 1 and the question of Pareto optimality (section 3.2.51. A 



summary of the demand properties is given in |3.3| 

We leave to the Appendix [A] the analysis of other kinds of pdfs, where we study the demand 



phase diagrams of IWP distributions with compact support (section A.ll and with fat tails 



(section A.2|. The demand phase diagram for a pdf with an arbitrary number of maxima is 



studied in A.3| — in details for a smooth multimodal pdf in Section |A.3.1| and on a simple 



example of singular bimodal distribution in section |A.3.2] — . 

Finally we summarize the main results and give several perspectives for further research in 
Section g] 
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2 Model of discrete choices with heterogeneous agents 
and positive externalities 

2.1 Agents model 

We consider a population of N agents (i = 1,2,..., N). Each individual i has to make a binary 
choice. Depending on the context, this binary decision may represent the fact of buying or 
not a good, adopting or not a given standard, adopting or not some social behavior such as 
joining a riot 03], or a journal club [77], [78], etc. Formally each agent i must choose a strategy 
uji in the strategic set = {0,1} (u>i — 1 denotes to buy /adopt /join, uji = otherwise Q 
Hereafter, without loss of generality, we will refer to the simplest market situation where the 
agents are customers who must choose whether to buy or not a single good at a price P. Our 
main concern is with the agents' behaviors, and P is considered as an exogenous parameter 
— e.g. it is posted by a monopolist selling the good — . Non-market models like those recently 
considered by, i.e., Glaeser et al. [30] are obtained by setting P = or by considering P as an 
exogenous social cost, common to all the agents. We are interested in the collective outcome of 
the agents decisions. In a following paper |44j focusing on the market context we will analyze 
the consequences of the customers collective behavior on a monopolist's program for fixing the 
optimal price. 

The population is heterogeneous. Each individual i has an idiosyncratic preference or will- 
ingness to pay/adopt (hereafter IWP) Hi, meaning that in the absence of social influences, an 
agent i adopts the state oji — 1 if Hi is larger than the price P. Following Mc-Fadden [62] and 
Manski [55], we work within the framework of Random Utility Models (RUM): we assume that 
the Hi are time independent random variables independently and identically distributed (i.i.d.) 
in the population. Denoting by H the mean and by a the variance of the IWP distribution, 
hereafter we assume that the random variable (Hi — H)/a is distributed according to: 

j^j jj 

V(x < — < x + dx) = f(x)dx, (1) 

a 

so that / is a pdf with zero mean and unitary variance. In the RUM view point, the agents 
have utilities H® and H\ for not adopting and adopting respectively, with H° and H 1 the 
corresponding population averages. For k = 0, 1, one writes H k — H k + e k , where e k are 
independent random variables with zero mean and standard deviation a k , with not necessarily 
identical pdfs for k — and k = 1. Then in our model we have Hi — H\ — Hf, H = H 1 — H° 
and / is the pdf of the normalized difference Xi = (e\ — ef)/cr, with a 2 = (a ) 2 + (a 1 ) 2 because 
of the additivity of the variances of independent variables. As particular examples, if e k are 
Gaussian variables, then Xi is also Gaussian; if e k are uniformly distributed on, say the intervals 
[— a k , a k ], then Xi has a trapezoidal pdf, which becomes a triangular distribution if e.g. a° = 
(no uncertainty in the utility of not adopting) . In the following we do not assume any specific 



form of the pdf of x. In Section 2.4 below, we present in details the class of pdfs considered in 
this paper. 

If all the individuals had the same IWP, the outcome in absence of social interactions would 
be very simple: either the price is below this common value, and everybody buys, or it is above it 
and nobody buys. All the individuals would behave in the same way, obtaining the same payoffs, 
and in the market aggregate analysis, they may be replaced by a fictitious representative agent 
[54] , In the case of a heterogeneous population considered here, only the agents with Hi > P 
would buy at price P, but getting different payoffs. 

The situation is more complex when the decision of each agent depends also on the decisions 
of others ( [59 and references therein) . We assume that each agent is the more willing to pay 



'Some authors use the notation s; = 1 and Sj = —1; both encodings are equivalent: it suffices to replace 
= (si + l)/2 in our model and identify the coefficients of corresponding expressions. 
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the larger the number of buyers in the population. We consider a linear separable surplus, that 
is if agent i buys at the posted price P, his surplus is 



Si = Hi + J V -P, (2) 
where rj is the fraction of buyers in the population. Taking into account the definition of ojf. 

N 

X 



1 N 



1=1 



We assume also that the externality Jrj corresponds to strategic complementarities, i.e. that 
the strength of the social influence is positive: J > 
The actual surplus of agent i is: 

Wi = Si on. (4) 

In order to maximize his surplus, agent i should buy/adopt (uji = 1) if Si > 0, but not (u>i = 0) 
when Si < 0. Since the IWP are i.i.d., when N is very large (more precisely, in the limit 
N — > oo), by the law of large numbers, the fraction of buyers ^ — which is the average of u>i — 
converges to the expected value of u>i over the IWP distribution. Thus, rj is given by the fixed 
point equation: 

ri = V(H t -P + Jr) > 0). (5) 

The marginal customer m, indifferent between adopting or not, is defined by the condition of 
zero surplus, S m = 0: 

H m -P + J V = (6) 

so that ([5]) may be written as 

n = V{Hi>H m ). (7) 
For what follows it is more useful to write ([5| as 

r) = V(Hi -H >-S) (8) 

where 

S = S{J, H,P- r]) = H - P+ Jr/ (9) 

is the population average of the (ex ante) surplus Si. It depends on the parameters J and H, 
which arc properties of the customers population, and P, the exogenous price. 

Notation. Generally, upon manipulating functions, we put in parenthesis the parameters 
in front, separated with a semicolon (;) from what we consider the variable. Our notation 
S(J, H, P; rf) in ^ indicates that J, H and P are considered as parameters, whereas r\ is the 
variable. Sometimes, when the context is clear, we drop down the parameters and keep only 
the variable, writing thus S(ry). Whenever we consider functions of two variables, they are 



separated by a simple colon, like in equations (37 1 and (38 1, and in Appendix A. 3 



2 More generally, the social term may be proportional to the fraction of buyers in an individual-depending 
subset of the population, called "neighbors" of agent i. In this paper, we consider a global neighborhood, where 
every agent has social connections with every other agent, mainly because this case can be studied analytically. 
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2.2 Aggregate behavior: normalized equations 

Clearly, the fraction of buyers rj depends on the strength of the social influence J, the price 
P and the average willingness to pay in the population H, and on the distribution of the 
deviation of the IWP Hi from its population average H . The agents choices depend only on 
the surplus sign, and they are invariant under changes of the surplus scale. Since the surplus is 
linear, we can formally multiply every term of the surplus by a same strictly positive number 
without changing the agents' choices. An adequate scale is given by the typical scale of the 
IWP distribution: it is convenient to measure each quantity (J, H, P) in units of the width a of 
the IWP pdf. Hence instead of four parameters, we are left with three independent parameters. 
Hereafter we will thus work with the following normalized variables 

J , H P 

j = - h=—, P =- 10 
a a a 

In addition, as it is obvious from equations ^ and ([9]), 77 depends on the price P and the 
average willingness to pay H only through their difference H — P. We introduce the normalized 
difference: 

H - P 

6= = h-p, (11) 

(T 

which is the average ex-ante surplus in the absence of externality. For short hereafter we call 6 
the bare surplus. In non-market models (p = 0) it is the average willingness to adopt. 



Remark. In (almost) all the following we will work with the above reduced variables (10 1, 



(111, referring to them as the (normalized) strength of social influence, average willingness to 
pay, price, and bare surplus. However one should keep in mind, especially when interpreting the 
results, that they represent the ratios of the non normalized parameters to the width of the IWP 
distribution. Clearly other normalizations are possible. An alternative of particular interest is 
the normalization obtained by measuring every quantity in units of the social strength J: the 
relevant parameters are then 

(tee-, h=—, P=-j, = — - — =h-p (12) 

(equivalently one can do as if J = 1). Note that this choice of normalization is no more than an 
equivalent representation of the parameters space; indeed one has a = h= h/j, p = p/j. 

It is also interesting to analyze the results in term of the set of parameters ( Jl2| ) , which allow to 
stress the model properties as a function of the degree of heterogeneity (relative to the strength 
of the social influence) — a homogeneous population corresponding to the limiting case a = 0, 
a highly heterogeneous one to a large a — . 



With the normalized variables (10 1, (111, equation d8| becomes 



r/= / f(x)dx = l-F(-s), (13) 

where F is the cumulative probability distribution and s = S/a, with S defined by ([9]), depends 
on h and p through the bare surplus <5, that is 

s = s(j,S;rj) =5 + jr). (14) 

If the pdf has infinite support, 

F(-s) ee V{x <-s)= f f{x)dx. (15) 



10 



In the case of a compact support [x m , xm], one can write: 

/>max{xjv/ . — s} 

V =l-F(-s)= / f(x)dx. (16) 

J max{i m , — s} 

Obviously, when — s < x m , we have rj = 1, and when — s > ir^f we have 77 = 0. 
2.3 Hints from two extreme cases 

In the absence of social influence the problem is simple because after introduction of j = 



in (14 1 we obtain s — 5 which does not depend on 77. Then, due to the monotonicity of 



cumulative distributions, the fraction of buyers ( 13 1 is a monotonically increasing function of 5 



(equivalently, at fixed h, a decreasing function of the price p): 

r) = l-F(-5). (17) 

Another extreme case is that of a homogeneous population: Hi = H for every i — a sit- 
uation obtained in the singular limit a = 1/j — > (the IWP distribution becoming a Dirac 
distribution) — . In that case every agent is faced to exactly the same decision problem, so that 
at equilibrium either 77 = or 77 = 1. For each agent the surplus in case of adoption would be 
H — P if no other agent adopt (77 = 0), whereas if 77 = 1 the surplus is H — P + J. In fact, 7/ = 
is a solution for H < P, while 7/ = 1 is a solution for H > P — J. Hence there is a domain, 
P — J < H < P. where the two solutions coexist. The whole population behaves as a single 
agent who either does not adopt, 7/ = 0, or adopts, r\ = 1, with a different surplus depending on 
whether he is "in" or "out of the market: this is analogous to the problem of multi-equilibria 
with hysteresis in trade analyzed by Baldwin and Krugman [B], except that here the problem 
arises only at the collective level. 

We have thus on one side, for J = and a finite, a unique well behaved equilibrium, and for 
J > but a = a situation of multiequilibria. The question addressed in the following aims 
at understanding what happens "in between" . We will show that when the social interaction 
strength is large enough compared to the heterogeneity width, the demand faces a complex 
problem. More precisely, there is a critical value jg of j — J jo. Below it, the fraction of buyers 



(equation (13)) follows monotonically the price variations. Beyond jb, equation (13 1 presents 
multiple solutions. Among them, the (possibly multiple) Nash equilibria are those solutions 
that have an economic meaning, i.e. for which the demand decreases when prices increase. 
The Section [3] of the paper is devoted to a detailed study of the nature of the solutions 



of equations (1131 and ( 14 1 with j > for distributions satisfying very general smoothness 



hypotheses, detailed in the next section. 



2.4 The idiosyncratic willingness-to-pay distribution 

Since we are interested in the generic properties of the model, we explicit the general charac- 
teristics of the idiosyncratic willingness-to-pay (IWP) distributions covered by our analysis. 

Since a pdf must be integrable, f(x) (equation ([lj ) must vanish in the limits x — + ±00. 
For sufficiently regular pdfs, this can happen in two different ways: either the pdf decreases 
continuously to as x — > ±00, or it is strictly zero outside some compact support [x m ,XM]- 
Most of the analysis in this paper is restricted to the class of pdfs obeying to the following 
hypotheses: 

HI. Modality: f is unimodal, that is it has a unique maximum 

f B =sapf{x). (18) 
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H2. Smoothness: f is non zero, continuous, and at least piecewise twice continuously diffcr- 
cntiable inside its support, ].r,„,arji/[ , where x m and Xm may be finite or equal to ±00. 
In the latter case / is stricly monotonically decreasing towards zero as x — > ±00. 

H3. Boundedness: the maximum of /, fs (that may be reached at x m or xm if these numbers 
are finite) , is finite: 

f B < 00 (19) 

Within the class of pdfs satisfying HI, H2 and H3, we will consider more specifically the 
important following prototypical cases: 

1. Unbounded supports: The support of the distribution is the real axis; the pdf is continuous 
and twice continuously derivable on ] — oo,oo[, with a unique maximum. A typical ex- 
ample, relevant to economics (see e.g. Q]), is given by the logistic distribution, although 
more generally we do not assume that the pdf is symmetric. We make the following 
supplementary hypothesis, that amount to impose that the pdf decreases fast enough for 
x — > ±00: 

H4. Mean value: the pdf has a finite mean value. Then, the smoothness condition H2 
imposes that / decreases when x — ► ±00 faster than | | 1 . 

H5. Variance: the pdf has a finite variance. Then, the smoothness condition H2 imposes 
that / decreases when x — > ±00 faster than \x\~ 2 . 

2. Compact supports: the support of the distribution is some interval [i mi xjn], with x m and 
Xm finite; the pdf is continuous on [x m , xm] and continuously derivable on ]x m , xm[, with 
a unique maximum on [x m ,XM\- Note that, since / has zero mean, x m < < Xm- 

Hypothesis H2 and H3 exclude cases where the pdf is not a function but a distribution — 
containing, e.g., a Dirac delta — . Clearly, if the pdfs support is the real line, ] — oo,+oo[, 
the boundedness hypothesis H3 is a consequence of the smoothness hypothesis H2. In the 
case of compact supports, H3 excludes pdfs diverging at a boundary of the support. Although 
hypothesis H3 is actually true under H2 if / is continuous on the closed interval [x m ,XM]> we 
explicit it because some of our results are valid under H3 even for pdfs less regular than those 
satisfying H2. 

Although hypothesis H5 is not necessary for the study of the aggregate demand, it corre- 
sponds to a wide family of realistic distributions for which one can conveniently use the standard 
deviation as the unit for measuring the relevant parameters (i.e., using normalization ( |10[ )). 

Generic results for unbounded support pdfs satisfying HI to H5 are presented in the main 
body of the paper. They extend previous results obtained for a logistic distribution [66 . The 
analysis of other types of pdfs is left to Appendix [XJ 

- In Appendix |A.1| we present general results for bounded support pdfs. The case of a 
uniform distribution on a finite interval [x m ,XM], which corresponds to an interesting 
degenerate case (/ is maximal at every point within the interval), has been presented 
elsewhere [45j . The particular case of general a triangular pdf is explicitly worked out for 
illustration. 

- In Appendix |A.2| we extend the analysis to fat-tail distributions, which correspond to 
an important limiting case of pdfs with infinite variance (for such distributions, the nor- 
malization constant cr in (jTj) and in equations (10 1 and (111 is no longer the standard 



deviation, but an arbitrary positive constant setting the units of H, J, P and C). 
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Finally, in Appendix A. 3 we extend the discussion to the case of multimodal pdfs (distri- 
butions with an arbitrary number of maxima) : we derive the demand phase diagram for 
a generic smooth multimodal pdf and we discuss the case of a singular pdf using as an 
illustrative example a pdf with two Dirac peaks. 



3 Aggregate choices and coordination dilemma 



In this section we discuss the demand function, 
fraction of buyers (or adopters, in non market 
we have already seen, this means studying the 
5 = h — p, and how it depends on the externality 
of values of the parameters j and 5, the demand 
as Nash equilibria from a game-theoretic point of 
the general hypothesis described in the precedinj 



that is the relationship between price p and 
contexts) 77, expressed by equation (13 1. As 



relationship between 77 and the bare surplus 
parameter j. We show that, for a large range 
presents two equilibria which can be qualified 
view. This result is valid for any pdf satisfying 
2; section. 



3.1 The direct and inverse demand functions 

The expected demand r\ d at a given value of S is obtained as the implicit solution of ( 13 1 and 



(14 1. As we will see, the application 77 — > 1 — F(—s(rf)) may be a multiply valued function 
of 77; it is thus preferable to express 6, or p = h — 8, as a function of 77, and determine the 
inverse demand function p (if), that is, the price at which exactly TV 77 units of the good would 
be boughtj^] 

Under the hypothesis H2 the cumulative distribution F is a continuous and strictly mono- 
tonic function on ]x m ,XM[ and has a unique inflexion point. Hence it is invertiblc. Denoting T 
the inverse of 1 — F(—s), we have the following equivalence: 



77= 1 - F(-s) 



(20) 



with s defined by ( 14 1. For unbounded supports, 1X77) increases monotonically from — 00 to +00 
when 77 goes from to 1 (see figure [I] for an example, and Appendix [A| for other cases). In the 
case of a compact support [ 

Xm, x m]> -T(^) fakes the finite values, T(0) — jcmi and T(l) — x 7n 
for 77 respectively and 1. Note that neither we assume / to be symmetric nor to have its 
maximum at x = 0. 

Replacing s in the r.h.s. of (20 1 by its expression ([Til) yields 



ith 



5 = V(j; V ), 



v(j; v) = r(77) - jn. 



(21) 



(22) 



Interestingly, T>(j; 77) depends on the parameter j but not on h. Actually, in Section 3.2 and in 
Appendix |A.3 we will have to consider T>(j\ 77) as a function of the two variables, j and 77. In 
the present subsection however, we consider T>(j; 77) as a function of the single variable 77, with 
j as a (fixed) parameter - hence according to our convention on notations introduced in Section 
2.1 we keep in mind the dependency of T> on j by writing T>(j; 77), and derivatives of T>(j\ 77) 
with respect to 77 are denoted T>'(j;rf). 

Plots of V(j;n) against 77 for different values of j are presented on figure [2] for the logistic 
distribution. Solutions to (21 1 correspond to the intersections of these functions with horizontal 
lines at y = S. 



3 In non market models, where generally p = 0, results in this section give the aggregate choice Nr/ d (j,h) as 
a function of j and h. This is the relationship between the fraction of adopters, the average willingness to adopt 
of the population and the strength of the social interactions. 
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Figure 1: T(rj) and derivatives as a function oft] for the logistic pdf of unitary variance. Remark: 
all these functions diverge at r) = and at r) = 1. 




Figure 2: V{j;rf} as a function of rj in the case of a logistic distribution of the IWP, for three 
values of j: j = 1 < j B , j = j B = 2.20532 and j = 5 > j B - 
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The stable equilibrium values of the demand satisfy 



V'(j;r 1 )>0. (23) 

Thus, the intersections of y — S with T> when V < correspond to unstable equilibria and will 
be ignored, as explained below. 

Given an ex tern ality parameter j and a given value of S, the solution 77 = i] d (j ; S) of equation 



(21 1 satisfying ( 23 1 gives the expected demand Nr] d (j; S) at a price p = h — S. 



From the definition of S and ( 22 1 , the inverse demand function is thus 

p d (r,)=h-V(j;r,). (24) 

This function depends on both parameters h and j, and when necessary we will write p d {rf) — 
p d (h,j;rj). 

As we will see, the demand r) d (j;5) can be a multivalued function of S for some range of 
parameters. On the contrary, since the function T(rj) is a uniquely defined function of 77, so is 
T>(j; rj). This is the reason why, instead of considering Jl3| , we prefer to obtain the properties 



of the demand 77 (j; S) from the analysis of equations ( |2l[ ) and (22 1. 

Under assumption H2, T(rj) is at least piecewise three times continuously derivable on ]0, 1[; 
its derivative T' is continuous and strictly positive. In particular, at any point 77 in ]0,1[, we 
have: 

V'(j;r,) = T'(r,)-j. (25) 

In the case of a compact support, the above equation also holds for the right and left derivatives 
at, respectively, 77 — and 77 = 1. 
In terms of the pdf /, 

T'(r,) = , with a = r(ry). (26) 

Under HI, V has a unique absolute minimum (qualitatively there is a unique point where the 
curvature of T changes from convex to concave; if T is smooth, it has a unique inflexion point). 
Thus 

minr'(77) = >0. (27) 

1 IB 

This minimum is reached at some value 77 = 

7/b = argminr'(r7). (28) 

v 

If / is smooth enough at its maximum, then T"{rjB) = 0: t\b is the inflexion point of T. For 
symmetric pdfs, t\b — 1/2, but we do not restrict to this case. 



As a consequence of the properties of r(?7), we see from equation (251 that T>'(j\ 77) is strictly 
positive for j < jb, with 

j B = r'( VB ) = i. (29) 

SB 

The value j B separates two regions where the model presents qualitatively different behav- 
iors. When j < js, the function T>(j; rj) is strictly increasing from —00 to +00 as 77 goes from 



to 1. As a result it is invertible: for any S in ] — 00, +oo[, equation (21 ) has a unique solution 

n d (6). 



If j > 3b, there is a range of values of S for which (21 1 has several solutions. 

The existence of multiple solutions in the demand is thus a generic property of discrete choice 
models with heterogeneous agents and social interactions (externalities). This is true whatever 
the number of maxima of /, as shown in Section |A.3| of Appendix [X] Actually, the domain 
where there is a unique solution, that is < j < js = l//s, is very narrow if f B is large: a 
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Figure 3: Demand r/ d (j; 5) as function of 5 = h — p in the case of a logistic distribution of the 
IWP, for three values of j: j = 1 < js, j = Ob = 2.20532 and j = 5 > OB- Notice that the 
origin of the horizontal axis (6 = 0) corresponds to h — p. Remark that prices increase from 
right to left. Unstable solutions: the negative slope curve joining {8l{o)i to {Su(J),f]u(j)) 

for j — 5 is not shown. 



leptokurtic distribution will have in general a narrower domain of unicity of the demand than 
a platykurtic distribution of same variance. 



In our case of unimodal pdfs, equation (21 1 may have three solutions for j > js (see 
figure [2]). The intermediate solution, laying on a branch with T>'(j; 77) < — where 77 increases 
as 6 = h — p decreases — is sometimes called a critical mass point in the literature |77j : it 
corresponds to a demand that would increase for increasing prices. Hence, in a tatonnement 
dynamics, this corresponds to an unstable solution separating the basins of attraction of the 
two stable equilibria. The marginal case, j — jg, is a bifurcation point (hence the subscript B) 



where multiple solutions to (211 appear on increasing j. The stable equilibria of the demand 
that satisfy ( 23 1 are represented against 5 on figure [3] 



3.2 Demand phase diagram 

3.2.1 The demand multiple-solution region 

Let us consider more in details the behavior of the application 6 — > r] d (d) in the case of a 
smooth unimodal pdf on ] — 00, +00 [. Considerations specific to compact support pdfs are left 
to Section |A.1| of Appendix [5] 

The functions r(7y) corresponding to pdfs satisfying HI to H5 are at least three times con- 
tinuously derivable on ]0, 1[, and diverge towards —00 and +00 as 77 goes to and 1 respectively. 
We have already seen that for j < jb there is a unique solution, and rf goes from to 1 as 
5 = h — p goes smoothly from — 00 to +00. 

For j > j B , (21 ) has 3 solutions whenever Su < S < 5 L (see figure |2j, where 5l and 6u are 



the values of 5 that satisfy the equality (marginal stability condition) in equation ( |23[ ). That 
is, the boundaries of the region with multiple solutions are the values for which T>(j;rj) has a 
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horizontal slope (see figure [2| : 

V'(j;r))=0 (30) 

which is equivalent to 

^M = 0. (31) 



Considering the definition (22) of V, this means that on these boundaries T>(j,7](j)), as a 
function of j, is the Legendre transform of T(r)). Under our hypothesis HI, T' has a unique 
minimum, and necessarily tends towards +oo as 77 goes to either or 1; T is strictly convex on 
]t]b, 1[, and strictly concave on ]0, t)b[, hence the Legendre transform is well defined and unique 
on each one of these intervals: equation p0| for j > js has indeed two solutions 77l(j) and 
Vu(j), given by 

j = T'( VA ), A=U,L. (32) 

with 

VlU) <Vb < Vu(j)- (33) 



From the knowledge of rju{j) and ??l(j), using (21 1 one gets the marginal stability curves 
&u{j) an d ^l(j), that is, the extreme values of 5 bounding the region where multiple solutions 
exist: 

Sa(j) = VVwaU)) = r(r/ A (j)) -3VaU), A = U,L. (34) 



As already stated, for 5u(j) < 5 < 5l{]), equation (21) has three solutions. The curve 
ri d (j; 5) has two stable branches (see figure]^]): an upper oner]f}(j; S) with ^(j; 5) > rju(j) > tjb, 
and a lower one 777^ (j; 5) with 777^ (j; S) < ^l(j) < ijb', they are joined by a branch of unstable 
solutions — the above mentioned set of unstable equilibria (see figure k3| — . The upper branch 

d d 

exists for 5 > Su(j), the lower one for S < S^(j). At the end points ^%-\l,u = 00. In other 
words, solutions with large fractions of buyers, i.e. high-r? solutions, only exist for 5 > 5u(j) 
whereas low-fy solutions exist only if 5 < 8l{])- Since Sjj(j) < <5i(j), the system has multiple 
solutions for the demand i] d whenever Sjj(j) < S < o~l{])- 

For j = js, these marginal stability curves merge at a single (degenerate) point SlUb) = 
SuUb) = 5b with 

( 5 B = -r'(7 te )? 7B +r(7 7B ) (35) 
This defines the bifurcation point B in the (J, 6) plane, 

B = {jb,S b }. (36) 

One should note that ?7(7,l(j) and 8jj,L{j), which depend on j (and on the function T(.)), are 
independent of h and p. 

3.2.2 Generic properties 

In fact, the preceding analysis can be made more general because the main results may be 
obtained only based on the continuity and the convexity properties of T, without assuming any 
smoothness properties of the derivatives of /. Let us consider this alternative. 

First, whatever the smoothness properties of /, the demand rf must be a decreasing function 
of the price: the economically acceptable values of the equilibrium demand, rf £ [0, 1], have 



to increase when 6 increases (p decreases). Thus, among the solutions of (21 1, the equilibria 



lie on the branches where D (defined by equation ( 22 1 ) , is an increasing function of 77 (for 
diffcrentiable pdfs, this condition is given by equation (23 1 ). 

Next, let us analyze T>(j; 77) as a function of 77 (sec figure |2|. By continuity of the function 
T(?7), T>(j;r]) is a continuous function of 77 6 [0,1]. As 77 — > 0, V — > —00, and as 77 — > 1, 
V — * +00. Since T is concave on ]0,77b], on increasing 77 from within [0,773], T>(j,r]) has 
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Figure 4: Demand phase diagram on the plane (J — J/cr, 5 — (H — P)/a), for a smooth 
IWP distribution (here the logistic). In the shaded region the demand presents multiple Nash 
equilibria. Outside this region, the demand is a single valued function of j and S. 



a maximum, ^l(j), on this interval. is by definition the Legendre transform of T(rj) 

restricted to ]0, tjb}- For 77 > rjs, T is convex, and thus T>(j, rj) has a minimum Sjj(j) on [r)B, 1], 
which is the Legendre transform of T(r]) restricted to [t)b, 1[- Beyond this minimum, T>(j;rf) 
increases with rj. 

Now, for j < js, the maximum on [0,775] and the minimum on [775, 1] of T>{j\rf) are both 
reached at ijb, hence T>(j;rf) increases monotonically as a function of 77 G [0,1]. Therefore, 



the solutions 77 to equation (21 1 are unique monotonically increasing functions of S for each j. 



As a result, the inverse demand (24 1 is a uniquely defined continuously decreasing function of 

77 e [0,1]. 

For j > js, the maximum 8l{j) is reached at 77 = t)l{j) G]0, 17b[. Beyond this maximum, 
V(j\rj) decreases as 77 increases. The minimum 5u(j) is reached at 77 = r]u(j) G]j7b,1[: there 
is an intermediate interval }r]L{j),Vu(j)[ containing 775 where T>(j;rj) decreases with 77, from 
5^ to Sjj. No value of 77 within this interval can be a stable economic equilibrium. Hence, 
for 5 ranging between these extrema of T>(j; 77) the demand i] d (j; S) as a function of 5 has two 
branches, a lower one for 5 < Sl, with 77^ (j; S) < r)L(j) < Vb and and an upper one for i5 > Su, 
with 7$ (j; 6) > r)u(j) > r] B . 

In the case of a continuously differentiable function, the preceding results are recovered, 
since the Legendre transforms — the above mentioned minimum and maximum of T>(j;rj) for 



j > 2b — arc reached at the values of 77 solutions of (30 1. All this discussion based on convexity 
arguments can be extended to multimodal pdfs, that is to cases where the distributions of the 
IWPs have more than one maximum. This is done in Appendix |A.3| 

3.2.3 The phase diagram 

The results of the preceding section may be summarized on a customers phase diagram in the 
plane (j,S), where we represent the boundary of the multiple solutions region, as in figure 
[4] These boundaries are the functions <5a(j), (A = L,U) defined by equations ( 34 1 , which 
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Figure 5: Demand phase diagram in the plane (er = a/ J, 5 = (H — P)/J), for a smooth IWP 
distribution (here the logistic). Inside the dashed region the demand presents multiple Nash 
equilibria. Outside this region, the demand is a single valued function of a and S. 
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are the two branches of the Legendre transform of T(rj), one for r\ < t\b and the other for 
77 > rjB- Note that in term of prices, the extreme values 8l(J) and 5u{j) correspond to prices 
PlU, h) < pu(j, h) given by 

p A (j,h) = h-8 A (j), A = U,L. (37) 

By construction of the Legendre transforms, the branch 8 = Sjj(j) is concave, and the branch 
8 = 8^(j) is convex. In addition, under the smoothness hypothesis, along each branch of the 
multiple solutions region in the phase diagram: 

^a(j) dV(j,r]A(j)) , Jr m 

d - = jj = -w), ; A e{L,U}. (38) 



This property may be easily checked by deriving (34) with respect to j and making use of 



(32 1. This means that the tangents to the boundaries have a slope given by the value of rj 
that is marginally stable on the corresponding boundary (i.e. by the rj value of the solution 
which appears/disappears as one crosses the boundary). Consequently, in the phase diagram, 
the width along the (5-axis of the multiple solutions region increases with j as a result of the 



convexity properties of the functions S\(j), (A = L,U). This may also be seen from (33 1, since 
the slope of the L boundary (the one corresponding to tjl) is larger than that of the U boundary 
(defined through rju)- At the bifurcation point B 7 these two boundaries merge, and, according 
to p8| ), have a common slope — rjg- 

Referring back to figure [2] upon increasing 6 from —00, we have the following picture: if 
j < js, the fraction r\ increases smoothly from and reaches its upper value 1 for 8 — > 00. 
That is, to each value of the bare surplus 8, — or each value of the average willingness to adopt, 
in non-market situations — corresponds a unique fraction of buyers/ adopters. In the phase 
diagram, figure |4j these solutions lie on the white region. On the other hand, if j > jg, when 8 
reaches the value 8u(j), a second, high-ry solution appears besides the low-77 one. These solutions 
co-exist for Su(j) < 8 < 5l{])- The low-77 solution disappears when 8 increases beyond 8l{j), 
leaving only the high-77 solution. The parameter values for which there are multiple equilibria 
is the grey region of the phase diagram, figure [4] Notice that in this region, there exists a third 
solution that we neglected because it corresponds to the unstable situation where the demand 
would increase with the price (or decrease with the bare surplus). 

As mentioned in Section 12.21 it is also useful to consider the same results in terms of the 



parameters a — a / J and 8 = (H — P)/J (see (12). The phase diagram in the plane (a, 8) 
is shown on figure [5j For large heterogeneity (a j J larger than ctb = there is a single 

smooth solution. For weak enough heterogeneity (ffg < 1/j_b), there is a domain with multiple 
solutions. In the limit a — > 0, one recovers the simple results for a homogeneous population, as 



briefly discussed Section 2.3 



3.2.4 Vicinity of the bifurcation point 

The vicinity of the bifurcation point B in the phase diagram is of particular interest. Under the 
smoothness assumption H2, we can study analytically its properties. Let us consider j = j's + e 



with < e << 1. Expanding (32 1 about ]b, remembering that T"{rjB) = 0, one gets, to the 



lowest order in e, the singular behavior 



<u = VB±J f ^- ) e 1/2 , (39a) 



SlmU) = r(»te) - Vb 3b ~ VBt T -J Y^—j e 3/2 - (39b) 
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The above singular behaviors are typical examples of scaling properties which are universal: 
the same scaling is obtained for any smooth distribution. From studies in statistical physics one 
expects the exponents (e.g., here, 1/2 for the behavior of rf) to depend mainly on the structure of 
the network of interactions: the exponents would be different at the analogous critical point for 
the model with agents situated on the vertices of a d-dimensional square lattice and interacting 
only with their nearest neighbors. Typically the exponents depend on d up to some critical 
dimension d c , above which they become identical to the "mean-field" exponents, which are 
those obtained here with the global neighborhood. For the present model, other universal 
scaling properties have been obtained, in relation with the hysteresis effects [80], and these 
have been used in order to analyze empirical socio-economics data |65j . In a related work (with 
a generalization of the model to more than two choices) , Borghesi and Bouchaud [TB] analyze 
empirical data for which the social strength can be estimated, and is found to be close to the 
critical value (the analogous of jg). 

In [44], where we consider the supply side, it will be seen that the bifurcation point B = 
UbiSb} m the (j, S) plane gives a singular point {jg, h B = S B } in the (j, h) plane which plays 
an important role in the phase diagram associated to the optimal strategy for the monopolist. 

3.2.5 Pareto optimality and coordination 

Each one of the equilibria r] d (j; S) discussed in the preceding section is a Nash equilibrium for 
the customers, at a posted price p. In this section we show that, whenever multiple solutions 
exist, that is for j > jg, the solution with the largest 77 is Pareto optimal. This is the solution 
n d (j;6) that satisfies 7] d (j;S) > r]u{j)- 

Let us recall that if a customer i decides to buy, it is because his (normalized) surplus 

8i = 5 + jri + Xi, (40) 

is positive. When s, < he doesn't buy and his surplus vanishes. Thus, his actual surplus is 
= s.iUJi (see Section [2|. 

Consider now the two equilibria r] d (j;6) in the region 6jj(j) < 5 < 5r,(j) (see the curve 
T>(j; rf) for j — 5 on figure^. Let's denote by 7?f (j; S) the low-77 equilibrium (r/f (j; S) < ?7l(j)), 
and by 77^ (j; S) the high-77 equilibrium (rju(j) < Vutii <5))- I n either equilibrium, the agents who 
buy are those with Xi > — S — jf] {j\ S). Since r}^ (j; S) < ^(j; S), agents with Xi < — S — jrjfj{j) 
are not buyers in any of the equilibria whereas agents with Xi > —S — jrj^j) are buyers in 
both equilibria. Those with —S — jr/^ij) < < — S — jrj^j) are buyers only in the high- 
77 equilibrium, and their utility is thus larger (strictly positive instead of zero) in that case. 
Moreover, even those agents that would buy in both cases have a larger surplus if the realized 
equilibrium is the high-77 onc - Hence, in the high-77 equilibrium all these agents have a larger 
surplus than in the low-77 one. This situation with two possible Nash equilibria, where the 
strictly dominant one may be risk dominated, is reminiscent of coordination problems in game 
theory. For a detailed analysis of this analogy see Phan and Semeshenko [74 . The present 
analysis shows that coordination problems may arise in systems with heterogeneous agents 
whenever the externalities are strong enough. 

Whether a Nash equilibrium — and which one in the case of multiple equilibria — is actually 
realized depends on the rationality of the agents and the information they have access to. 
In the context of bounded rationality and of repeated choices, a natural hypothesis is that 
agents estimate what will be the fraction of adopters, and may base their estimate on previous 
observations. In this paper we will not discuss these issues, that we are currently analyzing. 
Some partial results (dynamics with myopic agents and with various reinforcement learning 
paradigms) are discussed elsewhere [45] [79] . 
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Figure 6: Inverse demand p — h (= —5) as a function of rj for different externality strength 
values j, illustrated on the case of a logistic IWP distribution. 



3.3 Summary of the generic customers' model 

To summarize this section, if the social influence is small enough to satisfy the condition j < js, 
at each value of the bare surplus 5 — h — p, which measures the gap between the population 



average willingness to pay and the price, there is a unique solution r] d (j;5) to equation (21 1 
This demand is a monotonic increasing function of 5. However, if the social influence is large 
enough (j > js), there is a range of values 6jj(j) < S < Shii) for which two different (stable) 
solutions exist, a high demand one (i] d (j;S) > f]u(j) > rjs) and a low demand one (r] d (j;6) < 
T)l(J) < Vb)- In this region, the customers are faced with a coordination problem. If <5 is 
modified dynamically within this range, the demand may jump abruptly between these two 
solutions, a situation analogous to so called first order phase transitions in physics. Outside 
the range [Su(j)> ^i(j)]> there is a single solution, like in the small j case. 

We showed that the threshold js, which corresponds to the onset of a bifurcation in the 
customers phase diagram, is determined by the maximum fs of the IWP pdf: jg = l/fs- 
Although most of the analysis has been done for smooth pdfs, we have shown that the generic 
behavior stems only from convexity properties of the function T(rj), that is from the fact that 
the pdf f(x) is strictly increasing for x < xb and strictly decreasing for x > xb, where Xb is 
the mode of /. 

Some specific properties which arise for distributions with compact support are discussed in 
appendix |A.l and the case of distributions with infinite variance is studied in appendix |A.2| 



In appendix A.3| we extend the results of this section to multimodal distributions. 



In market contexts (and in particular for the market analysis of |44j). it is useful to consider 
the inverse demand p d instead of S, as in standard approaches. In figure [6] we plot the values 
of S = p — h as & function of the demand r\ and the strength j of the social externality for the 
case of a logistic distribution. 
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4 Conclusion and perspectives 



The model of collective behavior considered in this paper, under the general hypothesis detailed 
in Section [2] may be declined in both non-economic and economic contexts. In the first case, 
one is interested in the fraction of adopters, which corresponds to studying the demand function 
for an exogenous price in the second case. 

Like in many other models in the recent literature, we consider optimizing agents making 
binary choices, with willingnesses to adopt that depend additively on an idiosyncratic part 
(IWP) and on the choices of other agents. The population is intrinsically heterogeneous: the 
IWPs are drawn from some distribution of mean H and variance a 2 . In contrast with other 
most studied models, here the idiosyncratic willingness-to-adopt heterogeneity is frozen: it 
does not result from (time varying) random shocks. In other words, each agent's choice is 
deterministic with a well known (to him) IWP, and we concentrate on the aggregate behaviors. 
We analyze the equilibrium properties (Nash equilibria) characterized by the emergence of a 
collective behavior resulting from the combined effect of externalities and heterogeneity. 

Our results, for global uniform interactions — a global social influence of uniform strength — 
in the limit of an infinite population — through the application of the central limit theorem — 
are summarized on phase diagrams. The axes of a phase diagram are the relevant model 
parameters. In the case of the demand, the parameters are: S, the bare surplus (that is, in 
the economic context, the gap between the average IWP and the posted price), and j, the 
social influence strength, both parameters being measured in units of the standard deviation 
a of the idiosyncratic term distribution (see 2.2). In this space of parameters, one draws 
the boundaries between regions ("phases") of qualitatively different collective behaviors. The 
boundaries, where "phase transitions" occur, are lines of non-analyticity (e.g. the demand is 
discontinuous on the boundary) . In our problem, the main feature characterizing a given region 
is the number of equilibria in this region. 

From a constructivist point of view, our model encompasses the classical downward sloping 
demand curve as a particular case. Indeed, one of the main results for the demand is that, for 
very general IWP distributions, there is a region in the phase diagram with multiple equilibria. 
More precisely, if the IWP distribution is mono-modal, there are two Nash equilibria for any j 
larger than a distribution-dependent value js- For smaller values of j, the (Marshallian) demand 
curves are, ceteris paribus, downward sloping (i.e. monotonically decreasing with increasing 
prices). For large externality strengths (j > js), when the population average willingness to 
pay is small enough, the demand becomes not-monotonic (as in Becker's example [3]). This is 
a very general property of the model with additive externalities, and does not depend on the 
particular statistical distribution of the idiosyncratic preferences. We also discuss (although 
more briefly) the results for multimodal distributions — for which there exist several regions 
with multiple equilibria, with possibly more than two equlibria for some of them — , and present 
detailed analysis of many illustrative examples. 

An important contribution of this paper is to exhibit the detailed properties of the bound- 
aries of the regions (in the parameters space) where multiple solutions exist. These properties 
are generic, in that they depend only on qualitative features of the IWP distribution. 

Future work may extend the results presented in this paper in several directions. First, the 
individual preferences may include a stochastic (noise) term like in [36, 33, 9il ITUl H31 [HI IBT1 H] . 
on top of the idiosyncratic one. Second, the present paper concentrates on the equilibrium 
properties - that may be considered as the "static" analysis of long term equilibria in the 
Marshallian tradition. Further studies should focus on the process that makes the system 
reach one or the other of the possible equilibria. A first study, implying revision of beliefs 
in an repeated choice setting, has already been published |79j : in the region with multiple 
equilibria, interesting complex dynamics occur with a large family of different equilibria being 
reached, depending on the particular learning rule used by the agents. Third, literature on 
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marketing and studies of social psychology shows that choices very often depend on imitation 
effects or social influence. For example, the existence of externalities in the Communication 
and Information Technologies (CIT) sector is well established, and may result in a multiplicity 
of equilibria |76j . This may arise in other sectors also. Yet, empirical and econometric studies 
allowing identification of the corresponding preferences distributions and the strength of the 
social influence are lacking [T5]. Fourth, the influence of social networks topologies deserves 
further attention. Ioannides [50 reported results for the "Thurnstone model, i.e. homogeneous 
IWP and stochastic (logistic) choices, mainly for tree-like and one-dimensional networks with 
nearest-neighbor interactions. It would be interesting to explore how the phase diagrams for 
the model considered here — heterogeneous IWPs and deterministic choices — are affected by 
short range interactions. Analytical and simulation studies on the RFIM [SHI [73] show that 
the heterogeneity introduces hysteretic effects in the dynamics, with interesting path-dependant 
properties (return-point memory effect [80]). A statistical method to calculate the return points 
exactly, starting from an arbitrary initial state, has been recently proposed [82] for the simple 
case of a one dimensional periodic network with nearest-neighbor interactions (called cyclic 
topology in |50j). The impact of such properties on both individual and collective economic 
behavior remains to be investigated. 

In a companion paper [44] . we will discuss within the market context the profit optimization 
by a monopolist fixing the price. Other extensions of the results of this second part involve the 
study of how, with repeated choices, the long term equilibria depend on the entangled dynamics 
where customers and monopolist learn from each other. Finally, at least two directions are 
worth to be explored: the case of an oligopolistic competition and the consequences of Coase 
conjecture in the case of choices with externalities involving a durable good - an issue already 
addressed in the literature [61 , but not yet in the regime where multiple equilibria exist. 

Acknowledgements 

This work is part of the project "ELICCIR" supported by the joint program "Complex Systems 
in Human and Social Sciences" of the French Ministry of Research and of the CNRS. M.B. G., 
J.-P. N. and D. P. are CNRS members. This work has been partly done while M.B. G. and V. 
S. were with the laboratory LEIBNIZ, CNRS-IMAG, Grenoble. 



24 



References 

[1] S. P. Anderson, A. de Palma, and J.-F. Thisse. Discrete Choice Theory of Product Differ- 
entiation. MIT Press, Cambridge MA, (1992). 

[2] W. B. Arthur. El farol. Amer. Econ. Review, 84:406-, (1994). 

[3] W. B. Arthur. Complexity and the economy. Science, 284:107-109, (1999). 

[4] W.B. Arthur, S.N. Durlauf, and D.A. Lane Eds. The Economy as an Evolving Complex 
System II. Santa Fe Institute Studies in the Sciences of Complexity, Volume XXVII, 
Addison-Wesley Pub.Co., Reading Ma., (1997). 

[5] R. Artie and C. Averous. The telephone system as a public good: Static and dynamic 
aspects. The Bell Journal of Economics and Management Science, 4 (1):89-100, (1973). 

[6] R. Baldwin and P. Krugman. Persistent trade effects of large exchange rate shocks. Quar- 
terly Journal of Economics, 419:635-654, (1989). 

[7] P. Ball. The physical modelling of human social systems. Corn-plexus, 1:190-206, (2003). 

[8] G. S. Becker. A theory of social interactions. The Journal of Political Economy, 82 
(6):1063-1093, (1974). 

[9] G. S. Becker. A note on restaurant pricing and other examples of social influences on price. 
The Journal of Political Economy, 99:1109-1116, (1991). 

[10] G. S. Becker and KM Murphy. Social Economics. Market Behavior in a Social Environ- 
ment. Cambrige Ma. The Belknap Press & Harvard University Press, (2000). 

[11] R. Benabou. Equity and efficiency in human capital investment: The local connection. 
Review of Economic Studies, 63:237-264, (1996). 

[12] B.D. Bernheim. A theory of conformity. Journal of Political Economy, 102:841-77, (1994). 

[13] L. E. Blumc. The statistical mechanics of strategic interaction. Games and Economic 
Behavior, 5:387-424, (1993). 

[14] L. E. Blumc. The statistical mechanics of best-response strategy revisions. Games and 
Economic Behavior, 11:111-145, (1995). 

[15] L. E. Blume and S. N. Durlauf. Identifying social interactions: a review. Working paper, 
(2005). 

[16] C. Borghesi and J.-P. Bouchaud. Of songs and men: a model for multiple choice with 
herding. Working paper http://fr.arxiv.org/abs/physics/0606224, 2006. 

[17] J.-P. Bouchaud and M. Potters. Theory of Financial Risk and Derivative Pricing. Cam- 
bridge University Press, 2nd Edition, 2003. 

[18] J-Ph. Bouchaud. Power-laws in economics and finance: some ideas from physics. Quanti- 
tative finance, 1:105, (2000). 

[19] W. A. Brock. Pathways to randomness in the economy: Emergent non-linerarity and chaos 
in economic and finance. Estudios Economicos, 8:1:3-55, (1993). 

[20] W. A. Brock and S. N. Durlauf. Discrete choice with social interactions. Review of Eco- 
nomic Studies, 68:235-260, (2001). 



25 



[21] W. A. Brock and S. N. Durlauf. Interactions-based models. In Heckman J. and Learner 
E., editors, Handbook of Economics, volume 5, Amsterdam, (2001). North-Holland. 

[22] J. I. Bulow, J. D. Geanalopolos, and P. D. Klcmpcrcr. Multimarket oligopoly: strategic 
substitutes and complements. Journal of Political Economy, 93/3:488-511, (1985). 

[23] R. W. Cooper. Coordination Games, Complementarities and Macroeconomics. Cambridge 
University Press, (1999). 

[24] J. Crane. The epidemic theory of ghettos and neighborhood effects of dropping out and 
teenage childbearing. American Journal of Sociology, 96:1226-1259, (1991). 

[25] N. Curien and M. Gensollcn. Lcs theories de la demande de raccordement telephonique. 
Revue Economique, 2 mars:203-255, 1987. 

[26] G. Dcffuant, D. Ncau, F. Amblard, and G. Weisbuch. Mixing beliefs among interacting 
agents. Advances in Complex Systems, 3:87-98, (2001). 

[27] S. N. Durlauf. Multiple equilibria and persistence in aggregate fluctuations. American 
Economic Review, 81:70-74, (1993). 

[28] S. N. Durlauf. Nonergodic economic growth. Review of Economic Studies, 60:203:349-366, 
(1993). 

[29] S. N. Durlauf. Path dependance in aggegate output. Industrial and Corporate Change, 
1:149-172, (1994). 

[30] S. N. Durlauf. A theory of persistent income inequality. Journal of Economic Growth, 
1:349-366, (1996). 

[31] S. N. Durlauf. Statistical mechanics approaches to socioeconomic behavior. In B. Arthur, 
S. N. Durlauf, and D. Lane, editors, The Economy as an Evolving Complex System II. 
Santa Fe Institute Studies in the Sciences of Complexity, Volume XVII, Addison- Wesley 
Pub. Co, (1997). 

[32] S. N. Durlauf. How can statistical mechanics contribute to social science? Proceedings of 
the National Academy of Sciences, 96:10582-10584, (1999). 

[33] S. N. Durlauf. A framework for the study of individual behaviour and social interactions. 
Working paper, , (2001). 

[34] S. N. Durlauf. Groups, social influences, and inequality: A membership theory perspec- 
tive on poverty traps. In S. Bowles, S. N. Durlauf, and K. Hoff, editors, Poverty Traps. 
Princeton: Princeton University Press, (2006). 

[35] G. Ehrhardt, M. Marsili, and F. Vega-Redondo. Emergence and resilience of so- 
cial networks: a general theoretical framework. Working paper, ArXiv e-prints, 
http://arxiv.org/abs/physics/0504124, (2005). 

[36] H. Follmer. Random economies with many interacting agents. Journal of Mathematical 
Economics, 1:1:51-62, (1974). 

[37] S. Galam. A new multicritical point in anisotropic magnets. III. Fcrromagnets in both a 
random and a uniform longitudinal field. Journal of Physics C, 15:529-545, (1982). 

[38] S. Galam and A. Aharony. New multicritical points in anisotropic magnets, i. ferromagncts 
in a random longitudinal field. Journal of Physics C, 13:1065-1081, (1980). 



26 



[39] S. Galam, Y. Gcfcn, and Y. Shapir. Sociophysics: A mean behavior model for the process 
of strike. Mathematical Journal of Sociology, 9:1-13, (1982). 

[40] E. Glaeser and J. A. Scheinkman. Non-market interactions. In M. Dewatripont, L.P. 
Hansen, and S. Turnovsky, editors, Advances in Economics and Econometrics: Theory 
and Applications, Eight World Congress. Cambridge University Press, (2002). 

[41] E. L. Glaeser, H. D. Kallai, J. A. Scheinkman, and A. Shlcifcr. Growth in cities. Journal 
of Political Economy, 100:6:1126-1152, (1992). 

[42] E. L. Glaeser, B. Sacerdote, and J. A. Scheinkman. Crime and social interactions. Quarterly 
Journal of Economics, CXL507-548, (1996). 

[43] M. B. Gordon. An introduction to statistical mechanics. In Bourginc P. and Nadal J-P., 
editors, Cognitive Economics, pages 131-155. Springer, (2004). 

[44] M. B. Gordon, J. -P. Nadal, D. Phan, and V. Semeshenko. The perplex monopolist: Optimal 
pricing with customers under social influence. In preparation, (2007). 

[45] M. B. Gordon, J. -P. Nadal, D. Phan, and J. Vannimeuns. Seller's dilemma due to social 
interactions between customers. Physica A, 356, Issues 2-4:628-640, (2005). 

[46] M. Granovetter. Threshold models of collective behavior. American Journal of Sociology, 
83(6):1360-1380, (1978). 

[47] M. Granovetter and R. Soong. Threshold models of interpersonal effects in consumer 
demand. Journal of Economic Behavior & Organization, 7 (l):83-99, (1986). 

[48] W. Hildebrand. On the "laws of demand" . Econometrica, 51:997-1019, (1983). 

[49] J. Holyst, K. Kacperski, and F. Schweitzer. Social impact models of opinion dynamics. In 
D. Stauffer, editor, Annual Reviews of Computational Physics IX, pages 253-273. World 
Scientific, Singapore, (2001). 

[50] Y.M. Ioannides. Topologies of social interactions. Economic Theory, 28:559-584, (2006). 

[51] Ernst Ising. Beitrag zur Theorie des Ferromagnetismus. Dissertation, (1924). 

[52] M. L. Katz and C. Shapiro. Systems competition and network effects. The Journal of 
Economic Perspectives, 8(2):93-115, (1994). 

[53] R. Kindcrmann and J. L. Snell. Random Markov Fields and their Applications. American 
Mathematical Society, Providence, Rhode Island, (1980). 

[54] A. P. Kirman. Whom or what does the representative individual represent? Journal of 
Economic Perspective, 6(2):117-136, (1992). 

[55] B. Krauth. Simulation-based estimation of peer effects. Journal of Econometrics, forth- 
coming, (2006). 

[56] B. Krauth. Social interactions in small groups. Canadian Journal of Economics, 39 (2) 
May:414-433, (2006). 

[57] H. Leibenstein. Bandwagon, snob, and veblen effects in the theory of consumers' demand. 
Quarterly Journal of Economics, 64:2:183-207, (1950). 

[58] R. D. Luce. Individual Choice Behavior. Wiley, New- York, (1959). 



27 



[59] C.F. Manski. The structure of random utility models. Theory and Decision, 8:229-254, 
(1977). 

[60] C.F. Manski. Economic analysis of social interactions. Journal of Economic Perspectives, 
14 (3) Summcr:115-136, (2000). 

[61] R. Manson. Network externalities and the coase conjecture. European Economic Review, 
44(10):1981-1992, (2000). 

[62] D.L. McFadden. Conditional logit analysis of qualitative choice analysis. In Zarembka, 
editor, Frontiers of Econometrics, pages 105-142. New York: Academic Press, (1974). 

[63] D.L. McFadden. Quantal choice analysis: A survey. Annals of Economic and Social 
Measurement, 5 No. 4:363-390, (1976). 

[64] R. D. McKelvey and T. R. Palfrey. Quantal response equilibria for normal games. Games 
and Economic Behavior, 10:6-38, (1995). 

[65] Q. Michard and J. -P. Bouchaud. Theory of collective opinion shifts: from smooth trends 
to abrupt swings. The European Physical Journal B - Condensed Matter and Complex 
Systems, 47:151-159, (2005). 

[66] J.-P. Nadal, D. Phan, M. B. Gordon, and J. Vannimcnus. Multiple equilibria in a monopoly 
market with heterogeneous agents and externalities. Quantitative Finance, 5(6):557-568, 
(2006). Presented at the 8th Annual Workshop on Economics with Heterogeneous Inter- 
acting Agents (WEHIA 2003). 

[67] J.-P. Nadal, G. Weisbuch, O. Chcncvcz, and A. Kirman. A formal approach to market 
organisation: Choice functions, mean field approximation and maximum entropy principle. 
In J. Lesourne and A. Orlean, editors, Advances in Self-Organization and Evolutionary 
Economics, pages 149-159. Economica, London, (1998). 

[68] L. Onsager. Crystal statistics, I. A two-dimensional model with an order disorder transition. 
Physical Review, 65:117-149, (1944). 

[69] A. Orlean. Le role des influences interpersonnelles dans la determination des cours bour- 
siers. Revue Economique, vol. 41, n. 5:839-868, (1990). 

[70] A. Orlean. Bayesian interactions and collective dynamics of opinion: Herd behaviour and 
mimetic contagion. Journal of Economic Behavior and Organization, 28:257-274, (1995). 

[71] E. Ostrom. Collective action and the behaviour of social norms. Journal of Economic 
Perspectives, 14:137-158, (2000). 

[72] D. Phan, M. B. Gordon, and J.-P. Nadal. Social interactions in economic theory: an insight 
from statistical mechanics. In Bourginc P. and Nadal J-P., editors, Cognitive Economics. 
Springer, (2004). 

[73] D. Phan, S. Pajot, and J-P. Nadal. The monopolist's market with discrete choices and 
network externality revisited: Small- worlds, phase transition and avalanches in an ace 
framework. In Ninth annual meeting of the Society of Computational Economics, Univer- 
sity of Washington, Seattle, USA, July 11-13, 2003, (2003). 

[74] D. Phan and V. Semcshcnko. Equilibria in models of binary choice with heterogeneous 
agents and social influence, submitted to European Journal of Economic and Social Sys- 
tems, , (2007). 



28 



[75] J. Rohlfs. Theory of interdependent demand for a communications service. Bell Journal 
of Economics and Management Science, Vol. 5, No. 1, Spring:16-37, (1974). 

[76] J. Rohlfs. Bandwagon effects in High- Technology Industries. Cambridge, Ma., MIT Press, 
(2001). 

[77] T. S. Schelling. Hockey helmets, concealed weapons, and daylight saving: A study of binary 
choices with externalities. The Journal of Conflict Resolution, t. XVII, N. 3, (1973). 

[78] T. S. Schelling. Micromotives and Macrobehavior. W.W. Norton and Co, N.LY., (1978). 

[79] V. Semeshenko, M. B. Gordon, J.-P. Nadal, and D. Phan. Choice under social influence: 
effects of learning behaviors on the collective dynamics. In Topol R. and Walliser B., 
editors, Contributions to economic analysis. Cognitive Economics. New Trends, pages 177- 
203. Elsevier, (2007). 

[80] J. P. Sethna, K. Dahmcn, S. Kartha, J. A. Krumhansl, B.W. Roberts, and J. D. Shore. 
Hysteresis and hierarchies: Dynamics of disorder-driven first-order phase transformations. 
Physical Review Letters, 70:3347-3350, (1993). 

[81] C. Shapiro and H. Varian. Information Rules: A Strategic Guide to the Network Economy. 
Harvard Business School Press, (1999). 

[82] P. Shukla. Exact solution of return hysteresis loops in a one-dimensional random-field ising 
model at zero temperature. Phys. Rev. E, 62 (4) October :4725-4729, (2000). 

[83] A.R. Soetevent and P. Koorcman. A discrete choice model with social interactions; with an 
application to high school teen behavior. Journal of Applied Econometrics, forthcoming, 
(2006). 

[84] S. Solomon, G. Weisbuch, L. de Arcangelis, N. Jan, and D. Stauffer. Social percolation 
models. Physica A, 277:239-247, (2000). 

[85] H.E. Stanley Introduction to phase transitions and critical phenomena. Oxford University 
Press, (1971). 

[86] G. Topa T. Conley. Socio-economic distance and spatial patterns in unemployment. Journal 
of Applied Econometrics, 17 (4):303-327, (2002). 

[87] L. L. Thurstone. Psychological analysis. American Journal of Psychology, 38:368-398, 
(1927). 

[88] G. Topa. Social interactions, local spillovers and unemployment. Review of Economic 
Studies, 68 (2):261-295, (2001). 

[89] T. W. Valente. Network Models of the Diffusion of Innovations. Hampton Press, Cresskill, 
NJ, (1995). 

[90] B. Von Rabenau and K. Stahl. Dynamic aspects of public goods: a further analysis of the 
telephone system. Bell Journal of Economics and Management Science, Vol. 5, n2: 65 1-669, 
1974. 

[91] D. J. Watts. A simple model of global cascades on random networks. PNAS, Vol. 99 no. 
9:5766-5771, (2002). 

[92] W. Wcidlich. Sociodynamics: A Systematic Approach to Mathematical Modelling in the 
Social Sciences. Harwood Academic Publishers, (2000). 



29 



[93] W. Weidlich and G. Haag. Concepts and models of a quantitative sociology. Berlin, Hci- 
dclberg/New York, Springer- Vcrlag, (1983). 

[94] G. Weisbuch, A. Kirman, and D. Herreiner. Market organisation and trading relationships. 
Working paper 1996, published in: The Economic Journal, Volume 110 Issue 463:411-462, 
(2000). 

[95] G. Weisbuch and D. Stauffcr. Adjustment and social choice. Physica A, 323:651-662, 
(2003). 

[96] J. Zabel Y.M. Ioannidcs. Neighborhood effects and housing demand. Journal of Applied 
Econometrics, 18:563-584, (2003). 



30 




Figure 7: Triangular pdf of unitary variance and a maximum at xb = — 1. Left: V{rf) and its 
first derivative. Right: T>(j;rj) for different values of j . The values ofrjL, 5l, r\y and 8jj are 
represented for the particular value j = 5. 



A Appendix: Demand for other distributions 

In this Appendix we extend our analysis to more general distributions. 



We first ( A.l I explicit the particularities introduced on the above generic results when the 
pdf has a bounded support. In Section |A.2| we relax hypothesis H5, and consider pdfs with 
unbounded support called fat tail distributions in the literature. Finally, we extend our results 
to multimodal distributions in Section IA.3I 



A.l Pdfs with compact support 

We consider here pdfs f{x) with compact supports: x S [x m ,XM] presenting a unique maxi- 
mum (which may be located at one boundary). Clearly, such pdfs have finite variances. The 
discussion follows the same lines as that of the generic smooth distributions, except that in 
addition one has to pay attention to the values of V and its derivatives at the boundaries n = 0, 
77 = 1. A simple uniform distribution, analyzed in 45J, is a particular case where the maximum 
of the pdf is degenerate. 

In this section, derivatives like r'(l) or r'(0), stand for the left and the right derivative of 
r at 77 = 1 and 77 = 0, respectively. Due to the fact that the pdf strictly vanishes beyond its 
support, if the price is very high with respect to h (small 8) there may be no buyers at all, 
and 77 = 0. On the contrary, if the price is very low (large 8) the market may or not saturate, 
i.e. 77 = 1, depending on the behaviour of the pdf in the vicinity of x m . We represent the 
lines delimiting the regions where these solutions exist on the phase diagram. It should be 
stressed that these lines only indicate saturation and non-existence of a market. Their nature 
is different from that of the boundaries 5l and 5jj. In the following we consider a triangular 
distribution with a maximum fg at xb to illustrate our general results. The figures in this 
section correspond to xb = — 1) i.e. a case where the maximum lies inside the support. The 
case x b — — Xjyi has been considered in the study of the learning dynamics |79j . 

In the case of compact supports {x„i,xm], r(ry) increases from T(0) = — Xm < 0, to T(l) 
> 0. Like in the generic case of unbounded supports, T' reaches a minimum at t]b, and 
there is a critical value js = ^/fB beyond which multiple solutions appear. Notice that if 
the maximum of the pdf lies at x m or at xm, Vb lies at one of the boundaries of the [0,1] 
interval. As already shown in Section |3.1[ if the pdf is symmetric — as is the case for the 
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uniform distribution — , rjB = 1/2. 

In our example of a triangular pdf with maximum at xb we have: 



2(x- 



if x m <x<x B , 



I (x M -L)(xJ-x B ) lf *b<x<xm, 

with /(x) = outside the support. At its maximum, fg = 2/{xm — x m ). The constraint of 
zero mean and unit variance imposes a relationship between the support boundaries and the 
value of x B : x m = (-Xg - V 3 ( 8 - x 2 B ))/2, x M = (~xg + \/3(8 - x 2 B ))/2. 

If Xb = xm the distribution increases monotonically from at x m , reaching its maximum 
at xm, and is defined only by the first equation in ( |A-1 1. Conversely, if xb = x m then f(x) 



decreases inside its support, and is defined by the second equation in ( A-l ). In both cases, f(x) 
presents a discontinuity at one boundary of its support. 

The function T and its first derivative for the triangular distribution are respectively 

r / n = | -xm + y/ {xm - x m )(x M - x B )r) if < r/ < i] B 

I -x m - y / (x M - x m )(x B - x m )(l - rj) if rj B < V < 1 



and 



if rjB < rj < \ 



= < y ^-X-^) r. _ : ( a - 3 ) 



2- v /(l~'?) 

where 773 = (xm — xb)/(xm — x m ) is the inflexion point of T. They are represented on figure|7] 
(left) for the particular value xb = — 1- Notice that T"(rf) is discontinuous at rjB, because the 
maximum of f(x) is a cusp. 



The corresponding inverse demand function V(j; rj) (equation (22 ) ) is represented on figure 
[7] (right). Due to the finite range of the compact support, there are two new particular values 
of 5: 6 = T>(j; 0) = T(0) = -x M , independent of j, and 5i(j) = V(j; 1) = r(l) - j = -x m -j: 
for 5 < S , 1) = (no-market) is a solution, while for 5 > <h(j) there is a solution rj = 1 (market 
saturation). These extreme values of rj may be reached upon finite values of S (i.e. finite prices 
and finite average IWP) only in the case of compact supports. 

For j <Jb, 2?(j; v) is strictly increasing and invertible on ]0, 1[: for any 5 in ] — xm, ~x m — j[ 



equation (21 1 has a unique solution rj d (S) 7^ {0,1}. One can easily check that j < js implies 
j < x m — x m , so that —xm < —x m — j and consequently T>(j; 0) < T>(j; 1). In the particular 
case where / is the uniform distribution, one has precisely jg — xm ~ x m . For the triangular 



pdf (A-l ), ( |29[ ) gives jg — (xm — x m )/2, the support's half-width. Thus, for S < S < Si(j) the 
fraction of buyers/ adopters is a monotonic increasing function of 5. For 5 < do, rj — 0, and for 
S > Si(j) the market saturates (77 = 1). 

For j > jg there are two stable solutions whenever Sjj(j) < S < <5l(j)- Due to the existence 
of the extreme solutions 77 = and 77 = 1, the analysis is more cumbersome than for infinite 
supports. If the maximum of the pdf lies inside the support, the solutions r]L,(j) and T)u(j) 



of equation (30 1 lie in ]0,1[ and Su(J) an d ^i(j) both satisfy V = 0. On increasing 6 from 
—00, there is no demand until S = min{(5o, 5jj(j)}- If <5o < 5u{j), when 5 increases beyond 6q 
the demand becomes finite and remains unique provided that 6 < 6 < Su(j). For S > Su(j) 
we enter the region of multiple solutions. On the other hand, if 5 > 5jj(j), the system steps 
directly from the no-demand solution into a region where a finite demand equilibrium with 
rj > rjB coexists with the no-demand one. Notice that the high- 77 solution may correspond 
to either a fraction of buyers strictly smaller than 1 (if <5l(j) < Si(j)) or to saturation (if 
SlU) > Si(j)). In the case of the triangular distribution it is straightforward to check that the 
multiple solutions region sets in at jg = (i M — x m )/2, Sg = xg/2. 
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5=h-p 




Figure 8: Triangular pdf of unitary variance and a maximum at Xb 
diagram. 



-1: customers phase 



If the pdf has its maximum at one of the boundaries of its support, either r/u or t]l coincide 
with rjB- More precisely, if xb = x m then r\u = tjb = 1, if xb = %Mj r lL = Vb = 0. 

Summarizing, the customers phase diagram for pdfs with compact supports have two sup- 
plementary lines with respect to that with unbounded supports. They indicate the boundary of 
the viability region (no market exists below this line) and the saturation boundary (above which 
all the customers are buyers). Figure [8] presents the customers' phase diagram for our example 
corresponding to the triangular pdf of unitary variance ( A- 1 1 , with a maximum at xb = — 1. 



A. 2 Pdfs with fat tails 

Fat-tail distributions are characterized by the fact that the pdf / has a slow decrease at large 
values of x, so that the variance is infinite - or even the mean is infinite. Equivalently this 
occurs if r(ry) diverges 'too fast' to — oo when rj — > 0. In the case of the logistic, T ~ log 77; for 
the Gaussian, T ~ — \J— 2 log 77; for a power law, r ~ — \. This suggest to consider the general 
following smooth behavior: 



as 17 ~^ 0: T(r)) K(-logr)) a - 



(A-4) 



with the constant K > 0, a > and b > (ab 0). The fat-tail case corresponds then to a = 
and b > 1 in the above equation (A-4). 



A particular example of a fat tail distribution is a pdf with a power law decrease, which for 
large x behaves like: 



(A-5) 



with /i > 0. Then, for small 77, T ~ — \ with b = l//i, so that b > 1 means fj, < 1. For /i < 1, 
not only the variance but also the mean value of the random variable x is infinite. 



33 




Figure 9: T function and derivatives corresponding to the function (A-6) 



For fat-tails distributions one has to look at finite size effects: it is no more possible to 
take directly the large N limit and make use of the central limit theorem: quantities like 
(1/N) G{xi) for any function G will be dominated by rare events, that is by the largest 
values encountered in the population of (large but finite) size N . There is, however, no difficulty 
in doing this analysis: the results are obtained by doing as if the pdf had a finite support, the 
upper bound ijj being given as an increasing function of N (for an introduction to statistics 
with fat tails, see e.g. [T7]). 

Since we considered compact supports in the preceding section, we concentrate here on the 
marginal case fi = b — 1, which can be analyzed as a limiting case of distributions with infinite 
support. 

For fi = 1, f(x) does not have a finite variance. Then, the value of a that defines the 
normalized variables ( 10 1 may be any (finite) measure of the width of the distribution, as for 
example, the value of x at which f(x) is equal to half its maximum. Let us discuss this marginal 
case on a simple example (see figure [9| : 



r = - 



corresponding to the cumulative function: 



r\ 1 — rj 



(A-6) 



sgn(z) 



1 

4' 



(A-7) 



The corresponding pdf is, 



VT 



(A- 



as represented on figure 10 Since this is a symmetric distribution, r/s = 1/2, and one finds 
(see equations (28 1, (29 1 and ( 35 ) ) that the point B in the customers phase diagram is jb = 
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Figure 10: Pdf corresponding to equation (A-8) 



T'{rj B ) — 8, 5b = —4. Notice that, like for any monomodal distribution (satisfying thus 
hypothesis HI), T(?7) is convex for r\ > tjb and concave for r\ < r/s, with as before r/g being the 
inflexion point. The supply function has thus the generic behavior described in Section [3] even 
for fat-tail distributions. 



A. 3 Aggregate demand for multimodal pdfs 
A. 3.1 Smooth pdfs: generic properties 

In this section we consider the behavior of the application S — ► rj d (5) in the case of a smooth 
multimodal pdf with support on ] — oo,+oo[. The discussion Section 3.2 based on convexity 



arguments, can be extended to describe the phase diagram for the aggregate demand in the 
multimodal case. 

The minimal hypotheses we consider are the following. 



HAO. The pdf f(x) is continuous with a finite number K > 2 of x-values, 



-CO < Xn < 



< ... < x l B < 00 , for which / has a (possibly local) bounded maximum, 
f(x k B ) <oo fc = !,...,#. 



(A-9) 



For simplicity we assume also that the pdf is not constant on any interval of finite length. 
Actually, the discussion can be easily extended to less regular pdfs (in particular piecewise 
continuous pdfs), and pdfs constant on some intervals, but to keep the discussion shorter 
will leave that to the interested reader (in the case of a monomodal pdf, see the discussion 
on compact supports, and for the bimodal case see also below, Section [A. 3. 2| the singular 
case of a distribution composed of two Diracs) . 

• HA1. When considering smoother functions, we will assume / to be twice continuously 
diffcrentiable, so that in particular it has a zero derivative at every maximum and every 
minimum. 
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Let us denote by Xq the location of the minimum between a^ +1 and x k B . We assume / > 
everywhere on its support except possibly at some minima, and f{x) goes to zero as x — + ±00. 
By convention we set x'q = +00 and x^ = —00 (and we may write /(xq) = f(x^) = 0). 

In the monomodal case, we have seen that the critical value of j for the appearance of several 
solutions is js = 1//(xb). Here we will see that the relevant critical values are 

& = WT) k = 1 - K (A " 10) 



and also 



jr = —r-irT k=l,...,K-l (A-ll) 
3c f{x k c) 

(and it will be convenient to define as well j c = jf$ = 00) . 

Consider now the inverse demand at a given value of j. The ensemble of equilibria 6 d (rj) for 
77 € [0, 1] is the subset of the ensemble of solutions of ([21]) for which S increases (p decreases) 



as 77 increases. As for the monomodal case we study the function of 77 defined by (211 for 
any given j, 6(rf) = T>(j;rj). By continuity of the function T, S(rj) is a continuous function of 
77 € [0,1]. As 77 — * 0, S — > —00, and as r) — ► 1, 5 — ► +00. Increasing 77 from 0, 5(rj) increases. 
Similarly, decreasing 77 from 77 = 1, 6(rj) decreases. Since f(x) is continuous, r(?7) is continuously 
differentiable, with r'(?7) = dT(rj)/dr] = 1/ f(x) at x = —Y{rj). Hence r'(77) has (local) minima 
at values of 77 given by 

and (local) maxima at values of 77 given by 

T '^c) = TTk\ k = l,...,K-l (A-13) 

For a smooth enough pdf, the 77^ and tjqS are inflexion points for Y. Note that for any k = 
1,...,K, t?*- 1 < t?| < rj c . 

The most important remark is that T(rj) is strictly concave on every interval \rjQ~ 1 , rj B [, k = 
1,...,K, and strictly convex on every interval ]rj B , rj c [, k = 1,...,K. Then as 77 varies on 
[tj c ~ 1 , 77^], the function 2?(j, 77) = T(rj) — jrj has, at some value 77^ (j), a maximum 5* (J) which 
is by definition the Legendre transform of T(r)) restricted to [t^ -1 , rj B ]. Similarly, on [77^, t)q], 
T>(j, 77) has, at some value 7j£y(j), a minimum 8y(j), the Legendre transform of T(rj) restricted 

tO [77|, JJ*]. 

Depending on the value of j compared to the values j B ,jc, these min and max may be 
reached either at a boundary of an interval, or in the interior. More precisely: 

3<3% V k L=V k u = V k B (A-14) 

3% < 3 < 3c, V k B <Vu < Vc (A" 15 ) 

3%<3<3c~\ Vc' 1 < V k L < V k B (A-16) 

3c < 3, = rih = V k c (A-17) 

(and 77^ increases from 77^ to T]q as j increases from j B to j c , whereas 77^ decreases from ?] B 
to 77^, _1 as j increases from j B to j£. _1 ). In the case of a continuously differentiable pdf, every 
Legendre transform 77^ (j), A = L, U satisfies the marginal stability equation, 

One should note that 77^ L and L depend on j (and on the function T(.)), but not on h or p. 
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Now for j < 2b = mirifc j B , every min and max are reached at the corresponding value r/ B : 
this means that there is no intermediate regime with a decreasing behavior of 6(rj) as r; increases, 
hence S d (r]) — T>(j,r]), uniquely defined, is a continuously increasing function of i] <E [0, 1]. For 
3 > 3b, there is at least one k where the maximum 5^(j) is reached for rj = rj^(j) < T] B , and 
the minimum 5y(j) is reached for rj — Vi/U) 5* Vbi so that there is at least one finite interval of 
77 on which the function T>(j, rj) decreases with rj, and thus does not correspond to an economic 
equilibrium. Hence the demand rj d (6) has at least two branches. 

In the plane (j,S), the boundaries of the multiple solutions regions are thus given by the 
fonctions 6%(j) = T>(j,f]\(j)), A = L, U, which are the graphs of all the branches of the 
Legendre transform of T. By construction of the Legendre transform, every branch S — 6^(j) 
is a concave curve, and every branch S = S^(j) is a convex curve, and, under the smoothness 
hypothesis HAl, along each branch A = L, U, 



dj dj 



Recall that r)\ is the value of r\ for the solution which is marginally stable on this boundary. 

These boundaries can be easily drawn for any distribution making use of a parameterization 
by s (or equivalently x = — s): from the basic equations rj = l—F(—s) where F is the cumulative 
of the pdf /, s = T(r]), and T'(j]) = !//(— s); with the marginal stability condition ( |A-18 ) 



which gives j — T'(r]), the locus of marginal stability is then given in the plane (j,S) by the 
parameterized curve 

for x S support(/), 

j = l/Hx) (A-20) 

This is this representation that we have used to draw the phase diagram, figure |12| for the 
particular example of the bimodal distribution shown on figure |11| 

The domain of multiple solutions can then be described as follows. The phase diagram is a 
kind of superposition of diagrams associated to mono-modal phase diagrams, every maximum 
(every 'bump' in the pdf) k being responsible of the appearance of a domain of multistability: 
when j becomes larger than j^, a continuous solution split into two solutions, with a lower 
solution rj d (j,S) < uHj) < r\ k B and 5 < 5*1, and an upper one with r] d (j,S) > rj^(j) > r\ k B and 



<5 > <5^y (see figure 12 1. When j becomes larger than j^,, this bump is no more 'seen'. Since 
a minimum of the pdf, if not at a boundary, is in between two maxima, such an intermediate 
solution may exist either because of one bump or the other - or both. 

The branch 5 = (^j(j) has thus as left end point, B k = (j%,S B = T>(j B ,i] B )), and as right 
endpoint (if is finite), C k = (Jci"c = ^Uc'Vc))- ^ s ^ ne merging point of 5^ and S*, 
and Ck the merging point of 5^ and S^ 1 . Since 5 k and <5^ +1 must be both above 6y, these 
two branches must intersect one another for some value of j — j BC between j B and Jq. there 
is thus coexistence of three solutions in the triangular-like domain bounded below by 5^ (or 
max(i5^, <^ +1 ) if B k is below the branch (5^ +1 ), and above by S k for j < j BC , and by S k L +1 for 
3>j B c- 

In the smooth c ase (H Al), at every bifurcation point B k , resp. C k where two boundaries 
merge, according to A-19 there is a common slope —r] B , resp. —t/q . 



One may say that the pdf is probed at different scales for different values of j. Consider 
the graph y — fix). Every maximum below the line y — l/j is not seen (it does not change 
the structure of the solution), whereas a set of maxima higher than l/j, but joined by minima 
where / is still higher than l/j, is seen as a single global bump. This gives in particular that for 
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Figure 11: An example of bimodal pdf. 
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Figure 12: Phase diagram (aggregate demand) for the case of the smooth bimodal pdf shown on 
figure [77} 
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Figure 13: Examples of multimodal pdfs. At a given value of j = J/cr, the qualitative properties 
are obtained by looking at the intersection of the horizontal line y — l/j with the graph of the 
pdf y = f{x): for the particular value of j corresponding to the horizontal line on this figure, 
the two pdfs lead to the same qualitative properties of the Demand. 
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Figure 14: Phase diagram (aggregate demand) for the case of a bimodal pdf composed of two 
Dirac peaks. 



j > jg , the number of solutions is equal to one plus the number of times the line y = 1/ j cut 
the graph y = f(x) at points where / is increasing. Note that this does not give the number of 
solutions for a given value of S. On figure [13] two pdfs are shown; the intersection of the graph 
y = f(x) with the line y = 1/ j gives the structure of the demand at this particular value of j 
(in the case illustrated on the figure, the demand has 3 solutions for the two pdfs). 

A. 3. 2 A degenerate case: 2 Dirac 

Let us consider the particular case of an IWP distribution given by two Delta pics: Xi — ±x 
with equal probability (xq — 1/^/2 since the variance of / is normalized to 1). For j — 0, one 
has clearly r\ — 0, 1/2 or 1 depending on S < —x , ~xo < S < +x or S > xq. For j > 0, 
obviously r\ can still take only these three values. One gets easily the domain of existence and 



stability of these solutions, r\ = 0, 1/2, 1, by direct inspection of the equation ( 13 1 . The resulting 
phase diagram is shown on figure |14| 

This phase diagram for a singular distribution can also be understood by comparison with 
the predicted phase diagram for a continuous distribution. In the present case, the two maxima 
have equal height, +00, which gives j B = j B = 0, in agreement with the fact that boundary 
lines meet at j = 0. The minimum between the two maxima is at / = 0, hence jc = oo'- the 
domain of stability of the intermediate solution 77 — 1/2 extends to infinity, as it is the case 
whenever a minimum is at f(x^) = 0. The marginal stability lines are straight lines - hence, 
marginally concave and convex curves -, with slopes 0, 1/2 and 1 corresponding to the values of 



the solution marginally stable on the boundary, in agreement with (A-19). Since here there is 



no continuity in the demand at the singular points B 1 — (0, —1), B 2 — (0, 1), two branches do 
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not merge with a common slope: besides the fact that the demand can take only three values, 
this is the only place where the non smoothness of the pdf gives a feature of the phase diagram 
qualitatively different from what is obtained for a smooth pdf. 
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